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A COMPARISON OF ELECTRICAL AND DIFFUSION FORCES 
IN THE METABOLISM OF ELECTROLYTES 


ROBERT R. WILLIAMSON 
The University of Chicago 


Using an approximation method the electrical forces due to the 
metabolism of electrolytes have been compared to the ordinary dif- 
fusion forces and found negligible. 


In the development of mathematical biophysics it has been 
assumed for the purposes of simplification that electrical forces, cre- 
ated by a difference in the diffusion rates of positive and negative ions 
produced or consumed in metabolic reaction, would have a negligible 
effect on the concentration gradients of those ions. Preliminary con- 


siderations led N. Rashevsky to believe this to be true. By way of a 


preliminary study he considered the case of a spherical cell in which 


the diffusion coefficients are infinite, and only the permeabilities finite 


— 


(Rashevsky, 1938, p. 126, hereinafter referred to as M.B.). In an 
earlier paper (Rashevsky, 1935) he discussed the equations for the 
case of the consumption of two monovalent ions of opposite charge, 
present in equal quantities at an infinite distance from the cell in 
which they are consumed. Because of the difficulties in solving these 
equations the problem has awaited quantitative treatment by a new 
approximation method (M.B. Appendix). 

In this paper we shall apply an approximation method to the 
case of a spherical cell metabolizing to produce monovalent ions of 
opposite charge, which ions have different molecular weights and in- 
ternal diffusion coefficients but infinite permeabilities. Later we shall 
consider the cases where the permeabilities and external diffusion 
coefficients are finite and the cell is non-spherical. , 

In reasonably dilute solutions, where no forces other than those 
produced by random molecular collisions are considered, Fick’s law 
may be held to apply to the diffusion of metabolites. When, however, 


a field of force acting on the particles is present we must use a more 
- general formulation of the diffusion problem. 2 
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Let 

D = diffusion coefficient of an ion in cm? sec", 

c = concentration of an ion in gm cm”, 

B = mobility of an ion in sec gm", 

f = force on an ion in dynes, 

k = Boltzmann constant = 1-4 X 10-'* erg. degree", 
T = absolute temperature, 

J = rate of transport of an ion through the membrane 

in gm cm-? sec, 

q = rate of production of an ion in gm cm sec}, 
M = molecular weight of an ion, 

N = Avogadro’s number, 

é = charge on an electron in E.S.U, 

E = net charge on the cell in E.S.U, 

K = dielectric constant of the internal medium, 

V = volume of the cell in cm’, 

7’, = radius of the cell in cm, 


and let the subscripts + and — denote the ions to which the constant 
applies. Let a bar indicate “average,” and the subscript 0 “external” 
when applied toc. 

Then we may write immediately the following relations: 


D 

Bor? o 
=—Dgradc-+ Bfe ; (2) 

drt @- Cor Co- man 
Eee CosTe Sie f.=—f_. (3) 
If we consider that ¢ = c approximately at r,/2, then approximately 
—D grade = 4+ 2D-—, (4) 

0 


The total charges of the positive and negative ions considered sepa- 
rately are 


C. 
M. 


so that we may write the net charge on the cell as 


ad C. 
NeVand=-NeV : 


4 Soe 
B=sare(—) Ne. (5) 


The force on a positive ion at 7, is 
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Then from equations (2), (4), and (6), 


Co== Con An, N 2B; os C_ 


iP =" 2K ( M. va Cy ) 
(7) 
ap 2 Bypass tr 
Yo 3K M, M-_ 
The total transports and productions are 
Anrid,, AnrezJ.- and 4/327, q., 4/327° q-. (8) 


From considerations of a material balance in the stationary state, and 
from equations (7) and (8) we have, after rearranging 


Ne Di SN Bd) 
KkT M, * KkT M_ 


An?r,? 


Aanr —_¢2 Rina Genie, 
KET MS: KkT M. 
+ 6D_¢.— 1,2 ¢g-— 6D_¢-=90. 
Ne on 
Let 4 a7, KET p and, assuming K = 80 and r% = 10 em., we have 


p = 4.735 - 10%. Then dividing the first of these equations by CEP 
D,/M-_ and the second by ¢- p D_/M. gives 


(9) 

ae M,—- 6M. 6M. ena 
aie gates 0 

M_ p p C- 6D- 


Introducing (9a) into (9b) and solving for ¢, we obtain after rear- 
ranging 


+ Co- ). (db) 
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Le |= (91 10°) /6D+ + Cor __ 


P 
iM Be Oa iias A 
om “ED. to |% 
(10) 
9M_ 4.7 
255 - (a + Co+) C+ 
Oils Toe 
+ BF eu)? 


Since the linear and constant terms in ¢. contain the factor 10-” but 
in all other respects are similar to the square term coefficient, they 
will be negligible with respect to it. Then we have, after dropping 
them, 


Soe Es) OEY 1a Gy tee 6M, | 
Ca oM_ Cen + Co-) + oD) ( 6D. + Co+) te p ? 
or since 
1 q+ 10° May ORT 6" 6M, 
Ar tre gp es Pe 
we have 
C. =e M, ( q- To? qs. rT? Cox Co- 


1oM_D. | jeM.D. + eM. + oe - 


By means of equation (3) this may be written as 


4 = M, 
| ‘TD. eM. tT 
or simply 
Bait “ae Q+ To D,-+ D-_ i 
Resor ea a 


and similarly for c_ 


e ae Q- 7° D,+ D_ 
Te. ee (12) 


Here we may note the similarity of these equations to those of the dif- 
fusion of one substance as derived by the approximation method 
(M. B. Appendix). When D, >> D_, c, is proportional to ¢_ and to 
1/D_, and when D. >> D., ¢. is proportional to ¢ and to 1/D,. 
When (11) and (12) are introduced into (5) we get 


» ea 


PAN Ves 
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Papeol Ps ee Guts aD. -2D2Ts 
3 “or oe oes Oe 12) Mx 
Cae Co 
Poe |Ne=E, 
and introducing (3) 
B= 0, 


and the force vanishes with vanishing charge. 
If in equation (10) we let 


— a = c,? coefficient , 
”? 


O = €,. é 


c = constant term , 


Ln, 
a 
then we may put: 
£=2,+A, 
where A is a small quantity. Then «* — aa’ = —(ba +c). 


In taking « = a we have assumed —(ba + c) to be negligible. If we 
differentiate x? — ax? atx =a = x, we get (3%,.? — 2405) A a@, and 
to obtain a solution of equation (10) more exact than that given by 


(11) we must set this increment equal to —(bx + c) = — ab — bA 
ee. (oto — 202,) A — (3a? — 2a?) A= @ A= —ab—bA—C. 
ab + ¢ 
cee ae Piai 
For 
Meares Me =1074)5 
Deol On; DAS 
q.= 10° ; q-= 107 ; 
Hy =F peta 10> 
Cop 102-5 Coen LOTS: 
- and 
@=— 1.00, 10° 
P Eg el Be 
Cd DLs} 


then A, = —4.- 10-*° approximately. 
For the negative ions 
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0-—= 1:35) 2107 
= 3rd 0, 
Cae ees a One ee 


and A_ = 1.- 10-** approximately. 
E is then equal to 


or about 4.8 - 10+ #.S.U. on the whole cell. This is sufficient to ac- 
count for the observed charges on a cell (M.B., p. 126) yet it can be 
seen that it does not affect the diffusion gradients within the limits 
of experimental error. 

We may conclude therefore, that in a spherical cell the accumu- 
lation of a charge on a cell by selection of ions of different charge 
through different diffusion rates does not materially affect diffusion 
gradients or forces. We have also shown that an appreciable charge 
can be produced by metabolism using reasonable values of the con- 
stants of a cell. This may account for the charges that vanish with 
death of the cell as contrasted with charges produced by a Donnan 
equilibrium in the membrane, which would not vanish at death. 

In the special case just considered where h = D. = ow, the molar 
concentrations of the ions are identical at the surface, and since only 
average (linear) gradients are considered, any differences in molar 
concentrations will have the same sign throughout the cell; i.e., either 
c./M, > c_/M-_ for every r < 7 or ¢./M, < c_/M_. On the other hand, 
if D. and (or) h are finite, the molar concentrations need not be the 
same at the surface, and ¢./M. 2 c_/M_, depending on 7. This could 
produce a change in sign of the local charges as different regions are 
considered. The average charge may then be very small, but large 
local charges may exist which might create local forces of importance 
(in nuclear division, mitosis, protoplasmic streaming, or nerve excita- 
tion and conduction). Since these effects are of a purely local nature 
and since in the work that has been done previously (Young, 1939; 
Rashevsky, 1939) only average gradients, concentrations, and forces 
have been considered, the possibility of strong localized effects does 
not in any way affect the work of the past. The gross features of 
cellular mechanics will remain unchanged, and such considerations of 
local electrical forces need be made only when considering phenomena 
of a localized nature. The possible application of these local forces to 
such phenomena as nerve excitation and conduction will be discussed 
in the future. 


The author should like to acknowledge his indebtedness to Dr. N. 
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Rashevsky for suggesting this problem and for frequent helpful sug- 
gestions that have contributed materially to the paper, and to Dr. A. 
S, Householder for checking and correcting the manuscript. 

This investigation has been aided in part by a grant from the 
Rockefeller Foundation to the University of Chicago. 
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A CONTRIBUTION TO THE MATHEMATICAL BIOPHYSICS OF 
PSYCHOPHYSICAL DISCRIMINATION II 


H. D. LANDAHL 
The University of Chicago 


: The properties of a neural mechanism previously discussed are 
given in more general form. It is shown that, under certain condi- 
tions, a psychophysical scale value is proportional to a difference be- 
tween excitation factors at a synapse. More general relations are 
derived which hold among different variables in the two and three 
categories of judgment. An attempt is made to apply the results to 
the method of choices. Certain relations derived are compared with 
available experimental data. 


I 


In a previous paper (Landahl, 1938) the properties of a neural 
mechanism for discrimination were developed in analytic form by the 
use of an arbitrary distribution function of the fluctuation of excita- 
tion at a synapse. Let us now discuss certain properties without spec- 
ifying the distribution function. By so doing, it may be possible, to 
some extent, to correlate the results with some phases of psycho- 
physics. The results are then immediately extended. If any new re- 
lations are suggested in the process, it is so much the more fortunate. 
We confine our interest, however, to those situations in which it 
is possible, at least in principle, to obtain restrictive conditions upon 
the basic nerve parameters (Rashevsky, 1938, Ch. XXII) from” ex- 
perimental data. 

If (e —-j)’ is the value of the next excitation at s,* (Landahl, 
1938, Figure 1) due to random factors, we then have p (e—J)’ d(e—j)’ 
as the probability that the random fluctuation, (e—j)’, has a net exci- 
tation (stimulus) value which lies between (e—j)’ and (e—j)’ + 
d(e—j)’. As in the previous development [Landahl, 1938, equation 
(9)] we obtain for the probability of a wrong response, P,, when 
the two stimuli are presented, and only a correct or a wrong response 
_is allowed, 

* The following applies, irrespective of the symmetry of the function p, if 
we transfer the fluctuation to synapse 8,.- The transfer eliminates the use of 
P(—«x) [equation (2)], which if symmetry is assumed, equals 1 — P(x). 


—159— 


160 MATHEMATICAL BIOPHYSICS 


Py = { pe’ d(e—j)' (1) 
which, for the distribution used previously [Landahl, 1938, equation 
(7) ] reduces to the former expression [Landahl, 1938, equation (10) ]. 
However, if we do not specify the function p, equation (1) together 
with the definition of equation (2), 


P(x) = | p(e—i)' d(e—3) (2) 
becomes iss 
Pe= [-pd(e—i)'= [ “pd(e—iy'— [“pd(e—iy’ 
me e1-€2 -00 ~00 
Po 1 —P le, — 2.) (3) 
The probability for a correct response is then 
P, = P(é: — &2)). (4) 


Using the relations I, II, III [Landahl, 1938, p. 115], we may 
write for P’,, the probability of a wrong response, P,', the prob- 
ability of an equality response and P’,, the probability of a correct 
response, (the primes referring to the case in which the intermediate 
category is allowed) the following relations: 


Pro = | “pd(e—j)' =1—P(q—e, +h), (5) 
(€3-E2+h) 
‘ (€4-E2th) } 
Pa= | pdle—)' = Pl, —ee +h) —P(q—e—h), (6) 
(€3-E2-h) 
’ (€3-E2-h) : 
Pe= [ pd(e—3)' = P(e, —e%—h) . (7) 
The value of «, — e,, for the particular E(S) relation and 


p(e— 7)’ relation previously used, is given by [ Landahl, 1938, equa- 
tion (5) ] 


S Si — Sz 


sre t adel es lee (8) 
where 
AK 
L= one (1 — e*) (9) 


and A,K and a are nerve constants, ¢ is the time and S, , and S, are 
the stimulus magnitudes, S, > S,. The second approximate relation 
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of (8) holds when S, and S, are fairly nearly equal and is obtained by 
taking first order terms in the expansions. According to the approxi- 
mation one would be justified in plotting the P’s against per cent dif- 
ferences in the stimuli compared. 

If we now let one of the stimuli be the standard and permit the 
other to be larger than, equal to, or smaller than the standard, and 
let the subscripts g, d, and 1 mean “greater”, “doubtful” (“equal”), 
and “smaller”, we have P’, = P’a and 


Pie Piglet Pooh =P fOr oe sy 
S being compared with the standard S,, and 
Bisa P ys idiot oh Py Ep = Be L005 aie 


Then introducing 
S S—Ss 


x = Liogs-=L—<— , (10) 
we may write the relations 
Pees Pr), (11) 
Pr PID) (12) 
Pp he te (en it), (13) 
P= P(«a+h) —P(«x«—h), (14) 
PES Penh), (15) 


where now —o <x < o. From equation (9), « is a function of the 
time t, and the time relations are the same as previously discussed 
(Landahl, 1938). 

In deriving the above relations it has been assumed that the 
mechanism was symmetric. Then if P(x) # 4 for « = 0, so that a 
constant error (Culler, 1926, pp. 56-62; Guilford, 1936, pp. 27-29) is 
present, it is due either to a displacement of the distribution function 
or the fact that it may not be symmetric. Now if we assume that the 
function p is symmetric with a maximum at « = 0, but that the mech- 
anism (Landahl, 1938, Figure 1) is not, we have for the difference 
&, — &, instead of equation (8), the expression 


8, — 2 = —_ 1 — e-at) log ae 
: Nor 
Se ; 
Ai Bs (1 — e*) log ——. (16) 
Az hoe : 


Introducing S, Ss, and z, equation (16) may be written . 
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As K, SY Se soles —agt (17) 
i. 2 | 198 Fr Ne Noo Ja SE Dat 
or A a 
SY : 1=— 5 P 
9 a Leg NO aa ae cone (18) 
where 


iS A, K, Q (1 ae Caen) 


ase ee ee ee (19) 
PS AE Kren (lee a) 
hy ‘ 
¢'5 = Le log ——.. (20) 
Kee. 


We may then have a shift in the point of subjective equality or in the 
constant error due to asymmetry in the function p and also due to the 
asymmetry in the neural mechanism, ho; = ho. ; and the latter shift 
may change with time. 

If the thresholds of synapses s, and s, are not equal [ or similarly 
for s, and s, (Landahl, 1938, p. 109)] the resulting bias is equivalent 
to a shift in the maximum of the distribution function by an amount 
equal to the difference in the thresholds. Such an asymmetry in the 
mechanism would be a most plausible cause of the displacement of 
the maximum of p from x = 0. This type of bias might be termed a 
response bias as against a bias due to a difference in threshold of a sen- 
sory apparatus, ho, ~ ho. Consider, for instance, the case in which 
the left eye views the left stimulus and the right eye views the right. 
Let the left-right response bias, as discussed above, be eliminated. A 
bias might then be present due to a difference in sensitivity of the 


F AK : 
eyes, No, ~ hoo , or to a difference in the parameters Yes of equation 


(16)] of the nerve tract fibers (Landahl, 1938, Fig. 1, I and II). A 
measurable change in the bias might be obtained as a result of a one- 
sided change in the illumination conditions, introduction of color, the 
application of drugs affecting the pupil, the neural elements of the 
retina and the optic tract, ete. The information so obtained should 
be compatible with the present concepts of the parameters involved. 

We are discussing here only the case in which the stimuli being 
compared are presented simultaneously. If the stimuli are not pre- 
sented simultaneously, on the basis of the mechanism under considera- 
tion, the second stimulus would be very much favored. If the stimuli 
are not maintained together sufficiently long, no discrimination is pos- 
sible. However, by the addition of another mechanism to be discussed 
in a subsequent paper these temporal effects are changed. 

Let us now make the results more general by introducing an un- 
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specified relation «(S) between «, the value of the excitation factor 
at a higher center (as s; or S,), and S, the external stimulus magni- 
tude. Following an argument similar to that above, we obtain-equa- 
tions identical to equations (11) to (15) but with x defined by 


x = e(S) —e(Ss) , (21) 


instead of by equation (10), S being the stimulus which is compared 
with the standard Ss. 

The probability P, of a response S greater than Ss is, according 
to equation (12), equal to P(x), where x is given by (21). The pro- 
portion of responses of S greater than Ss may be determined empiric- 
ally under given conditions and, if the conditions are suitable, we may 
assume that the proportion will equal the value of the probability, 
P(a). Then, since P(x) is known, if the distribution function, p, is 
known, or may be assumed to be Gaussian, « is determined in terms 
of some unknown unit. As long as the unit cannot be determined in 
terms of the synaptic threshold h,, (% and « have the same dimension 
as the threshold), it might well be expressed as a standard deviation 
unit according to the customary methods of scaling (Guilford, 1936, 
p. 216). Under these conditions, a scale value is proportional to a 
difference between two excitation factors at the discriminating cen- 
ter. 


at 

Let us now discuss the relations (11) to (15) in connection with 
available experimental data. We first shall discuss the data by F. M. 
Urban on lifted weights (Urban, 1908). Seven subjects were asked 
to judge seven weights with respect to identical standards 300 or 450 
times each and the judgments allowed were “heavier”, “guess heavi- 
er”, “guess lighter’, and ‘“Jighter.” The “equality” or “doubtful” cate- 
gory was taken as the “guess heavier” plus “guess lighter” judg- 
ments.* 

If we identify the various proportions experimentally derived 
with the respective probabilities, and if we assume that the probabil- 
ity distribution function is Gaussian, or 


ul {(e-j)’ — 201? 
Ce (22) 


p(e— J) = ee 


*In Urban’s notation the symbols used were respectively h, hg, lg, and 1, 
with T —hg + lg. Then in transferring to the notation used above we have 
iy SS ae hel gph Fg a aac P,,and T= Pq. The change in nota- 
tion nade so as to have a consistent notation for all the experiments discussed 
and in order not to conflict with the notation previously adopted. 
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then from equations (15), (12) and (14), P’,, P, and P’, + P'a are 
normal cumulative frequency curves. In order to show the agreement 
with experiment in a most sensitive manner the experimental points 
are plotted on an arithmetic probability scale, on which the points 
should then fall along three parallel straight lines. This holds since 
the differences S — Ss are small compared with Ss, and from equa- 
tion (10), x is very nearly linear with S — Ss. 

In Figure 1 the averaged results from seven subjects are plotted 
on arithmetic probability paper. The results indicate essentially three 
parallel straight lines. The individual plots are similar to the average 
in this respect. The centrally located data, P, are slightly displaced 
toward the right of an exact central position. If this is taken as sig- 
nificant here, as it definitely is for individual cases, the interpretation 
on the basis of the mechanism is that the rise in threshold h for the 
three category case is not the same at s; and s, (Landahl, 1938, Fig. 
1). Such a bias does not seem unnatural. The result is the same as if 
there were a shift in the mean of the distribution function from the 
two to the three category case. The standard deviation of the distribu- 
tion functions, « , was taken as 5.7 grams. The mean for the two cate- 
gory case was taken as 97.1 grams and for the three category case, 
97.4 grams. The standard stimuli weighed 100 grams. The value of h 
was taken as 2.1 grams. Here, as elsewhere, we take the value of x 
as if it were given in the physical units of intensity. Strictly, from 
equation (10) we should involve the coefficient L. However, the unit 
of x is still not known so, for simplicity, we refer to « as if measured 
in grams, etc. 

On examining the visual and auditory data (Kellogg, 1930), we 
find somewhat more complex results. Before discussing the data, let 
us more specifically examine the meaning of the two and three cate- 
gory cases on the basis of the mechanism under consideration. We 
shall assume that generally h ~ 0 in equations (13), (14) and (15) 
so that the three category case is the more general situation. Now 
when the two category case is required, h must necessarily be zero 
when the judgment is given and hence we shall assume that the stim- 
ulus “two category” or its equivalent leads, through some center, C, 
to an excitation at s; and Ss. (Landahl, 1938, Fig. 1) which reduces the 
threshold to zero. 

It might be that even in the three category case there may be a 
tendency for the stimuli to influence the center, C , which lowers the 
threshold of s; and s,. This would happen if excitatory fibers from 
s, and s, (Landahl, 1938, Fig. 1) lead to C. Then, if the para- 
meters a of the fibers leading to and from C are large compared with 
those of the other fibers of the mechanism, and if their thresholds are 
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neglible, the increase in excitation at s, and s,, or decrease in h, is a 
function of «, — 7, (Landahl, 1938, Fig. 1). Now «; — 7, is a monoto- 
nically increasing function of «, and « and in the steady state, the 
two quantities are proportional. With a proper choice of the nerve 
parameters, A/a = B/b = 1 (Rashevsky, 1938, Ch. XXII), we may 
make the constant of proportionality equal to unity. Then, since 
& — & = |x|, (6: > &, see Landahl, 1938), when S, and S, are re- 
placed by S and Ss [equations (8) and (10)], A is a function of |a|, 
h(|zx|). 

Since ¢, — & = || acts as a stimulus to the fibers leading to the 
center C, the amount of excitation ats, and s, due to C would in- 
crease directly as the expression 1 — e-7l"| [Rashevsky, 1938, Ch. 
XXII, equation (5)]. If hk, is the initial threshold and], 6, and a are 
constants, then 


h=h see (l—eelel) , (23) 


or if h tends toward zero for large |x| , then 
h = h, e-!"1 , (24) 


Now h(a) may be determined in the following manner. From 
equations (6) or (14) we have 


(w+h) 


Pa= | pd(e—y)’. 
(#-h) 
Differentiating with respect to x we have 
ad P'4 d d 
ay ee Aan tid i el Cha hg Died AC ome fips AS) (25) 
Expanding, we obtain 
dP’; 
Te = [Pe the @) +50") +] + | 
— |p) —hp'(2) +p" (2) +. |j1— ed 
(26) 
=2hp'(2) +2p(2) S 4 hep" (a) 42 
Integrating, we have 
Pa =2hp(x) + fh? p" (x) dh-+.. 
(27) 


SiGe ee ae eee 
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Solving for h(x) we have 


iP: he ‘ 
h(a) hopper LEAT OE AM) Acai (28) 


where the second term on the right is a correction term which gives 
an estimate of the error made in the approximation. Since P’; and 
(x) may be determined empirically, h(a) may be obtained. In the 
inset of Figure 1 the h(a) relation is shown for lifted weights, h and 
a being both given arbitrarily in units of grams. In this case, h and x 
are essentially independent. The correction indicated by the second 
term on the right hand side of equation (28) was made. The maxi- 
mum correction was 0.10 grams or 5%, and the mean absolute cor- 
rection was 0.05 grams or 2.7%. The next higher order corrections 
are much smaller and neglible. 

In order to determine p(x) in equation (28) the following meth- 
od was used instead of the more direct but less reliable method of 
measuring the slopes of the empirical P(x) relation. Since the P(x) 
relations are so nearly Gaussian, the assumption of normality may 
be made and hence, for a given x, p(x) is determined, as is also 
»" (x). Actually, the P, = P(x) value is here given equal weight with 
the average of the other two points for each 2, P’, and (1—P’). If 
this is done, a graphic average is necessary, as otherwise the result 
will be systematically increased when the P values are less than 50 
per cent, and decreased when greater. As an example, the numerical 
average between 1 and 9 per cent is 5 per cent as compared with the 
graphical average of 3.3 per cent. 

If h is not independent of x, plotting P’,, P,, and P’, + P’4 on 
arithmetic probability paper will not result in, three parallel lines. In 
Figures 2 and 3 are shown averaged results of visual and auditory 
data by W. N. Kellogg (1930) and in neither case can three parallel 
lines (or even curves with parallel tangents at each point) very well 
approximate the data. In both these sets of data, the two category and 
three category cases were presented separately. Each of five subjects 
made 288 judgments on two sets, each set including both the two and 
the three category method. In the averages given in the figures not 
all sets could be included as the intervals were not uniform. For the 
visual data represented, each point was based on 720 j udgments while 
for the auditory data each point was based on 1152 judgments. 

As may be seen in Figures 2 and 8 the h(x) relation is not con- 
stant. In the insets of the figures are plotted the relation as derived 
from the data by means of equation (28). The correction term was 
taken into account as in the case of the lifted weights. The order of 
magnitude of the correction term in per cent was the same in the 
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three cases. As in the lifted weight data, the results for individual 
cases in both the light and the sound experiments were distinctly like 
that of the respective averages shown in the figures. 

In order to fit the curves for the visual data (Kellogg, 1930) it 
was assumed that h(a) was given by the derived equation (24). The 
curve in the inset of Figure 2 was obtained from (24) by plotting the 
logarithm of h against « and determining the slope and ordinate in- 
tercept of a straight line drawn through the points, thus obtaining 
ho = 0.49 meter candles and a 6 = 1.14 meter candles*. In each 
of the cases x is measured from the mean, 0.10 meter candles. The 
three curves in Figure 2 were obtained as follows. The central straight 
line, (P,), represents a normal cumulative curve having a standard 
deviation, o«, of 0.58 meter candles and with the mean, z,, at 0.10 
meter candles the line should pass through the solid circles. The up- 
per and lower curves were obtained respectively from equations (15) 
and (13) with the appropriate values of h given by equation (24) 
and with the constants given above. The curves should pass through 
the open circles and crosses respectively. The x values were again 
measured from the mean. The standard intensity was 21.68 meter 
candles. 

The ordinates in Figure 2 are differences between the standard 
and the comparison stimuli in meter candles. We are justified in using 
such a scale since |S — Ss| << Sg and the approximate relation in 
(10) holds accurately. However, the particular ordinate scale used 
in this case is due to the experimental procedure employed. The sub- 
jects were asked to respond “left” or “right” according as they judged 
the left or right stimulus to be the darker. In our notation, P; repre- 
sents the proportion of judgments “left” and P, represents the pro- 
portion of judgments “right”. The sign of the ordinate is then rather 
arbitrary. In the figure a positive difference means that the stimulus 
to the left is the darker. The constant error was such’ that, for equal 
stimuli, the one to the left was judged the darker. 

In the derivation of equation (24) the threshold of the fibers ex- 
citing the center C was neglected. Had this not been done, the curve 
in the inset of Figure 2 would have had a flat top to a distance, on 
either side of the mean, corresponding to the threshold neglected. The 
curve would then continue as shown. However it should be recalled 
that the same randomness present at the discrimination center would 
also be superimposed on the h(a) relation derived in equations (23) 
and (24). The effect would be to round off any discontinuities other- 
wise present. The resulting trend is illustrated in the figures by brok- 
en lines (Figure 2). 

In the experiment on sound stimuli (Kellogg, 1980), the sub- 
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ject was asked to report “stronger” or “weaker” according as he 
judged the comparison stimulus to be stronger or weaker than the 
standard. Since in this case S — Ss may be almost one-fourth of the 
value of Ss the linear approximation in (10) cannot be used as the ap- 
proximation breaks down rapidly when the difference is much over 
ten per cent. When a linear ordinate is used, the graph corresponding 
to Figure 3 shows an average trend which is noticeably concave up- 
wards. The use of the more exact logarithmic ordinate as in Figure 
3, and in accordance with equation (10), makes it possible to obtain 
excellent agreement with the theoretically derived relations. The out- 
side scale is the natural logarithmic scale. The inner values are the 
sound intensities expressed directly in (millivolts)?. 

The h(x) relation in the case of the auditory data, plotted in the 
inset of Figure 3, is definitely asymmetrical. In order to fit the data 
it is necessary to allow for the asymmetry by introducing two values 
for the parameter a 6 of equation (24). For S > Ss orx > 0,08 is 
taken as 0.35 ordinate units while for S < Ss orx <0,a6 = 2.8. 
Also h, = 0.93 in terms of the ordinate units or approximately 5 (mil- 
livolts)?. The distribution is assumed normal and the standard devi- 
ation, o , is taken as 0.18 or approximately 9 (millivolts)?. The three 
curves in Figure 8 are then obtained from equations (15), (12) and 
(13) in which h is given by equation (24), with the parameters as 
given above. 

From Figure 3 it may be noted that there is some evidence that 
the standard deviation of the data from the two category method is 
smaller than that of the three category (Kellogg, 1930, p. 43). 

From a comparison of the h(x) relations plotted in the insets of 
Figures 1, 2 and 3, one is lead to conclude that the form of the rela- 
tions in each case is perhaps significantly distinct, especially since 
the individual cases resemble their respective averages. Let us sup- 
pose that they are in fact distinct and that this result is not due to a 
chance grouping of the 17 subjects involved.* The differences may 
then be due to the modalities involved or to the experimental proce- 
dures or both. From the brief description of the experimental proce- 
dures given above it is clear that the procedure for, the lifted weight 
experiment was distinctly different from that for the sound or light 
experiments. In an experiment on lifted weights by E. Culler (1926) 
with himself as subject and with the same procedure followed as that 
of F. M. Urban (1908), there is an indication that h might decrease 


* Four of these subjects could not be included in the averages under discus- 
sion for reasons explained above. Though the number judgments for the indi- 
vidual cases is not too large, 144 to 450, the trend of the h(«) relation empiri- 
cally derived is very stable, especially in the visual and lifted weight data. 
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with |x|. But in the experiment by J. P. Guilford (1936) on one sub- 
ject, with constant errors experimentally eliminated and with the 
two and the three category cases separately presented, there is rather 
distinct indication of an h(x) relation similar to that of the visual 
data. The experimental points are in this case rather scattered. There 
is a-suggestion, however, that the mode of presentation is of impor- 
tance in this particular respect. If this suggestion should prove to 
be correct, it will require a natural explanation on the basis of the 
mechanism. 

For each of the five subjects in the visual experiment the values 
of oc , ho and oa 6 of equations (22) and (24) were obtained in the units 
used above. Between h, and a 6 a correlation of 0.1 was found. The 
correlation between hy) and o was 0.9, and between a @ and o, 0.2. 
These correlations and others, though not reliable, do impose restric- 
tions on the mechanism under discussion and they must be satisfied. 
For the present we shall not attempt to interpret these not too reli- 
able results. 

In Figure 2 it is clear that the three category data do not lie 
along two straight lines. The question then arises as to the meaning 
of the standard deviation unless, as above, a variable h is introduced 
If the mean slope is taken, the value obtained will vary with the 
range. Evidently a numerical determination of the standard deviation 
is no longer satisfactory since the initial assumption of normality is 
not too good. This is particularly the case in determining the critical 
ratio of a difference. Two curves which are clearly different may have 
a neglible critical ratio, and conversely, two curves having the same 
mean slope (on probability paper) may show a fairly large critical 
ratio obtained numerically. 

In this section the results obtained in the previous section are 
compared with experimental data. Certain revisions in the initial 
assumptions were necessary in order to obtain agreement with the 
data. The purpose has been to illustrate the results obtained, to dem- 
onstrate the method of determining the parameters and to interpret 
on the basis of the mechanism the meaning of any significant lack of 
agreement with experimental data. 


III 
Let us now attempt to apply results previously obtained to the 
Wahlmethode of Fechner. Consider again a structure similar to that 
of Figure 1 (Landahl, 1938) but having » first order fibers instead 
of two, and each first order fiber exciting n—1 inhibitory fibers so 
that there are n(n—1) inhibitory fibers instead of 2, Let n stimuli 
be presented simultaneously and let the stimulus magnitudes S; be 
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unknown but be so ordered that S; > S; fori <j. If Pi; is the prob- 
ability that a response is made to S; as the largest stimulus when all 
nm stimuli are presented and if P.i; is the corresponding probability 
corrected for chance, then 


n 1 
perro Gs tae, . (29) 


Pos = 


From equation (46) of the previous development (Landahl, 
1938), we have in the above notation 


Lé iS Sa" 
n—1 Bom erga: 


where L is defined by equation (9) and & corresponds with the recip- 
procal of o? of equation (22). We shall assume & to be independent of 
n.* P.., is the probability corrected for chance that the largest stimu- 
lus S, is responded to as the largest. Then from equations (29) and 
(30) we may write 


Poa = (30) 


te pe 3 SS 
Pete Spee SORE ra (31) 
Now define z; by 
2,=Lé log 3 , (32) 
So 


where S, is an arbitrary value. Comparing equation (32) with (10) 
we see that x is proportional to z; for S = Si and Sy = S). Equation 
(31) then becomes 


1 Legs 
=—+2,—-— - 33 
Py 5 z. Pea (33) 
Now let S, be removed and the remaining n—1 stimuli be pre- 


‘sented. Then in place of equation (33) we may write 


Le) ee (34) 


Ree GES 


Poo = Herat 


or more generally 
Nee 


1 1 : 
Oy 7) ee ee é 35) 
ean aay n—it+l 2% (20) 


~ * Tn the previous paper (Landahl, 1988), compare £ in equation (46) with k 
_ (Landahl, 1938, p. 110) in equations (57), (58), and (18) and also with the ex- 
pansion of equation (19) together with (20) and (18), . 
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The probability P., that S, is responded to as the largest, when 
there are stimuli, is less than the probability P.. that S. is responded 
to when the largest stimulus S, is removed and S, is the largest of 
those left. Let us assume that the former, P;, , is equal to the latter, 
P.», reduced by the factor (1 — P,,), or that the relative values re- 


main unchanged, so that 


P2, = (1 — P,,) Po» (36) 
and more generally 
Pi; = (1— Pyj) Pier - (37) 


From equation (37) we have 
Pi, = (1— Pi) (1 — Pros) +++ (1 — Poi iny) Pi (38) 
which in turn may be written because of equation (37), 
Py = (1 — Py — Poy — +++ — Poin) Pix. (39) 


Substituting from equation (35), the above equation may be written, 
if we define 


3; 237 Sle. (40) 
k=1 k=7 


I 
Pi, = —P, Spa ae iced (DE \ Mipcom ot ae er 
1 (1 1 Py ED as i 1 
(41) 


x [a+ Se+ (n—i+ 1) a — 3. 


From equation (41) we have n—1 equations in n unknowns, 2; , 
since the n-th equation reduces to an identity 


1 
P= (1 — Py, — ++» — Poy4y1) I (1 + 31 — 2@n + 2n — DS). 
Therefore introducing 
<3 Dt 
Yin = %j — aa (42) 
or 
Dy, 
ij = 24 — ——— 
Vij reer TEE (43) 
so that 
Sve Oy, (44) 


we have for equation (41), 
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Yin = Pu 1 2S 1 
i (1 — Py, — +--+» — Puia:) Sas eRe ) (45) 


k=1 


or because of equation (39), we may also write 


1 4-1 
Via =Ps— (3 tm +1). (46) 

If the n stimuli are presented and the proportion of responses 
that S; is the largest is obtained for each 7, or what is equivalent, Pi, , 
then the values of y;; are given by successive substitutions in the equa- 
tions represented by (45). Or if choices are made with the n stimuli 
_ present giving among other values P,, , then the most frequently chos- 
en stimulus, S,, is removed so that the next group of choices deter- 
mines P»» , etc., the values yi: , may be most readily obtained from the 
equations represented by (46). For a given value of L é, the values 
of S; are given by equation (32) if the values of z; are known. Or if & 
is known and also the values of z;, we may obtain # values. From 
equation (42) the values of z; are seen to be the same as the values of 
Yi,, but displaced by an unknown amount 3,/n, a value which de- 
pends upon the stimuli used. The origin of measurement of 2; is then 
not determined without other information. 

It should be noted that the equations do not imply that Yin > 
ee if Ps > Pe, or the converse. The possibility of inversions 
arises because it is assumed that the relative distribution of responses 
with and without the largest stimulus is the same in the derivation of 
equation (37). A stimulus for which the corresponding value of Yi 
is inverted would then be expected to be most associated with the 
stimulus just larger or smaller than itself. The presence of a particu- 
lar stimulus, then, not only changes the 2; or Yi: values uniformly but 
may do so non-uniformly. 

We have assumed that ¢ is independent of n. Depending on the 
details of the mechanism other relations may hold. If é is found to be 
independent of n, the use of the above procedure would be justified. 
However, its use would be restricted to the purpose of obtaining an 
additional relation among the parameters and effect of stimuli upon 
each other. The procedure suggested by J. P. Guilford (1937) leads to 
essentially the same results as that of the first section, as far as the 
evaluation of parameters is concerned. 

On applying certain results previously obtained to the method of 
choices it is found that only one restriction can be imposed upon the 
parameters; that the scale values derived may not be correctly or- 
dered, and will generally be displaced by an unknown amount; and 
that the absolute unit will not be obtainable. 


176 MATHEMATICAL BIOPHYSICS 


The author is indebted to Dr. A. Weinberg for checking the cal- 


culations and reading the manuscript. 
The publication of this paper was made possible by a grant from 
the Rockefeller Foundation to the University of Chicago. 


REFERENCES 


Culler, E., 1926. “Studies in Psychometric Theory.” Psychol. Monog. 35, No. 163, 
pp. 55-1387. 

Guilford, J. P., 1986. Psychometric Methods. McGraw-Hill Book Company, New 
York and London. 

Guilford, J. P., 1987. “Scale Values Derived from the Method of Choices.” Psy- 
chometrika, 2, 139-150. 

Kellogg, W. N., 1930. “An Experimental Evaluation of Equality Judgments in 
Psychophysics.” Arch. Psychol. No. 112. 

Landahl, H. D., 1988. “A Contribution to the Mathematical Biophysics of Psy- 
chophysical Discrimination.” Psychometrika, 3, 107-125. 

Rashevsky, N., 1938. Mathematical Biophysics. Chicago: The University of Chi- 
cago Press. 

Urban, F. M., 1908. “The Application of Statistical Methods to the Problem of 
Psychophysics.” Philadelphia, Psychol. Clinic Press. 


BULLETIN OF : 
MATHEMATICAL BIOPHYSICS 
VOLUME 1, NUMBER 4 
DECEMBER, 1939 


STREAMING IN DIFFUSION FIELDS 


GALE YOUNG 
The University of Chicago 


Two simple systems are discussed in which diffusion of a dis- 
solved substance causes continued circulation of the solvent, thru 
the mediation of gravitational forces in one case, and of solute viscos- 
ity forces in the other. The first is analagous to the circulation in a 
liquid heated from below. 


Living cells frequently exhibit a steady streaming of their proto- 
plasm. This phenomena has been the subject of considerable observa- 
tion and experimentation, but there is as yet no satisfactory theory 
of the mechanism involved. Any such motion may be regarded as the 
result of driving forces operating against the viscous resistance of 
the protoplasmic substance, but the nature of the driving forces is 
not known. They appear to be intimately connected with the cell 
metabolism, and one suggestion is to look for their source in the pro- 
cess of diffusion, as has proven quite promising in studying the me- 
chanics of cell division. 

With these considerations in mind we shall discuss some cases of 
solvent motion in a diffusion field. 


I 

Consider a pair of vertical tubes of length L connected to form a 
loop and filled with a solution. Suppose the concentration of the dis- 
solved substance to be maintained constant at the top and at the bot- 
tom of the loop, the values of the concentration being c” and c’ re- 
spectively, so that there is continual diffusion along the columns of 
solution. For simplicity we shall regard the tubes as one dimensional 
systems and deal only with average densities, velocities, pressures, 
forces, etc., over a cross section. The distance x will be measured up- 
ward from the bottom, capital and small letters will refer respectively 
to solvent and solute, the velocities V and v will be positive upwards, 
and the two tubes will be distinguished by subscripts 1 and 2, number 
1 being chosen as the tube where the solvent is flowing upwards. | 

The gravitational pressure acting downward upon the solvent 
column in a tube is 


—177— 
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g | Caz, (1) 


g denoting the gravitational constant and C the concentration, while 
the diffusion drag acting upward is 


L 
Dp a) [cae (2) 

C0) 
where ” and »’ are the solute pressures at top and bottom. The total 
pressure tending to make the solvent circulate up on side No. 1 and 


down on side No. 2 is therefore 


¢,) dz, (3) 


°L 
J=g | (CAB gee 
0 


since p’ and p” are the same for both tubes and hence drop out in the 
subtraction. The quantity J is the driving force which acts against 
the viscous resistance of the solvent. 

The flow must be the same across any cross-section, hence we 
have the equations of continuity 


CY AG, 
(4) 
Cre be 
where K and k are constants. The solvent flow (unlike that of the 
solute) will be equal and opposite in the two tubes, 

K = K, =—K, ; (5) 
and the viscous resistance to flow in a tube is proportional to the rate 
of flow. Thus, equating the driving force to the resistance, 

J = OLK (6) 


where © embodies the solvent viscosity and the cross sections of the 
tubes. 

In general a solvent expands slightly when a solute is added and 
its density is correspondingly decreased. For dilute solutions this may 
be adequately represented by the linear relation 


C == Go = be , (7) 
where b is a constant and C, is the density of the pure solvent. Ex- 
pression (3) then reduces to 


JG Cree Bs [e—e) de. | (8) 
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It remains to find the solute distribution in the two columns. Ob- 
viously if the solvent is at rest the situations will be identical, we shall 
have c, = ¢,, and there will be no driving force. If, however, there 
is circulation of the solvent, then c, + ¢c, and J may be different from 
zero and such as to maintain the circulation. This is related to work 
of Griffiths (1898, 1915). 

Assuming (Young, 1938) a drag force between solvent and solute 
proportional to the product of their densities and their relative ve- 
locities, the solute pressure equation is 


a = —ge —fcC(v—V) ; (9) 


which, introducing (4) and (7), and the Van’t Hoff relation 


aC 
(10) 
ele: 
aM 
becomes 
dc 
—_— = : Ah 
= = A+ Be; (11) 
wherein 
4a PBC 
104 
(12) 
K Ci 
pa tenld +f g 
104 
This has the solution 
: A-+ Bc = Ge, (13) 


where G is the integration constant. The boundary conditions c = ¢ 
atx =0,c=c" atx = L give 


A+ Be =G 
(14) 
A-+ Be” = Ge ; 
from which 7 
¢c’ e —s Cc” 
A = B Terre ’ (15) 


this being essentially a relation between K and k expressing the fact 
that the rate of flow of the solute depends upon that of the solvent, 
as is apparent physically. From the above relations it follows that 
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- Be ___ @BL 
c—c" = (c’ —c’") = ’ (16) 
and ; = 
L ; e 
0 


The bracketed expression in (17) is a monotonic decreasing func- 
tion of BL , as may be seen by considering 


€ 
o(y) = i= 


We find ¢(—«) = 0, ¢(0) = —1. The derivative function 
dp yroe—(#—1) 
dy = y(e’—1)? 


does not vanish at y = 0, and any other root would require a non- 
zero solution of y/2 = + sinh (y/2) which is impossible. Hence 
d ¢/dy is everywhere negative. 

Expression (8) may be written as 


J=g(1—b) (c’— ec") L [o(B;, L) —$(B, L)). (18) 
Now tube No. 1 has been chosen as the one in which the solvent moves 
upward, so that K, = —K, is positive. The solute flow will be great- 


er on the side where it moves with the solvent, and, whichever the 
sign of ce’ — c”, we shall have k, > k.. Then from (12) it is seen 
that B, > B,, and hence ¢(B, L) < ¢(B,L). Hence J is opposite in 
sign to c’ — c” , and since J must be positive in (6) it follows that c” 
must be greater than c’ for the driving force to be in the right direc- 
tion to tend to maintain the motion. Thus the solute must be diffusing 
downward for steady solvent circulation to be maintained. The solu- 
tion is then less dense at the bottom, and the phenomena is analagous 
to circulation in a liquid heated from below, and the results and meth- 
ods of this latter problem are of service in the diffusion case (Ray- 
leigh, 1916; Jeffreys, 1926; Low, 1929; Hales, 1935; Schmidt, 1935, 
1937). The case of time variable diffusion along a moving column 
has been treated by Somers (1912). 

It is apparent that the system is in equilibrium when the solvent 
is. at rest, for then the solute distribution is identical in the two col- 
umns. It will be shown, however, that the system is unstable if the 
solvent is stationary and the concentration gradient c” —. ¢' is suf- 
ficiently large. 

In (12) kis a rather complicated function of K , but since the 50- 
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lute flow will be aided by increased convection in the’ same direction 
it is seen that /, increases monotonically with K,, while k, decreases 
monotonically. Hence as V increases from zero, B, increases and B, 
decreases so that J is monotonic with K. If dJ/dK at K = 0 is great- 
er than OL of (6), then any small circulation will tend to increase 
and the state of zero circulation is unstable. 

Negiecting the second term in the numerator of (12) we have 
dB,/dK = —dB,/dK = f/a, so that from (18) at kK = 0 


dJ Me es 
ee) ACO) (19) 


_where ¢, is the value of d $/dy at Yo corresponding to K = 0. The 
inequality sign enters because the term neglected in (12) would in- 
crease the value of dJ/dK. At x = 0 we find by series expansion that 


¢ = —1/12; 4,’ will be of this order of magnitude. Thus the state of 
zero solvent circulation will be unstable if 
6) 

(P25 LS x (20) 


279 (1 —b) ($0) * 

Estimating ¢.’ to be —1/12 this is approximately 
6a0 

fg(1—b)° 


The lowest stable circulation rate may be estimated by expand- 
ing ¢ by means of the Bernouli series, convergent for 0 < w< Ax, 


il 1 il ay Hi 
Sabha as eS etal oh ge CP rs Ra 292 
S15 et aT TB 30.24 7 ere 


(226) > (21) 


Expression (18) then becomes 
b , ” L? B B (B,? — B,*) L? 
J=9(1— ye ory i 1 — B,) + 


_or, neglecting the b term in (12), 


++ | (28) 


LT? L? 2 k2 2 | 
= 90d) —0') S24 S Fe +35) en 


Equating this to OLK and removing the root K = 0 gives ~ 
60027, g? L? 640 
He al 20a? g(1—b) Lf(e” —¢’) 


as the approximate equation for the lowest stable circulation rate. 
As c” — c’ increases this expression tends to a limiting value 


i? (25) 
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which, with g/a neglected, is 


) a 
K=+T7.1—-. 26 
ELE Fr. (26) 
Now f is proportional to T/D , where D is the ordinary Fick diffusion 
coefficient and 7 is the temperature, so that with a as given by (10) 
it is seen that (26) varies as D/ML, M being the molecular weight of 
the solute. 


ai 

We now consider a system in which the solvent driving force 
arises not from gravity but from the viscosity of the solute. Suppose 
the tubes of the circuit to contain some sort of porous structure rigid- 
ly attached to the walls and thru which the solvent and solute have to 
move, such as a grid or network or a sponge structure. This gives a 
backward drag opposing the solute motion which will be assumed 
proportional to the total rate of solute flow, namely @ Lk per column 
where @ includes the solute viscosity and a structural factor. Instead 
of (3) and (8) we now have 


SS Op ks eke). (27) 
while (9) changes to 
dp 
ep ele =—0CDe. (28) 
This is the same as (11) except that now 
k 
a ee (6 a FCe)e 
(29) 
= i (K+ bk). 


As before, k is a rather complicated function of K, but k, in- 
creases monotonically with K, while k, decreases. Hence if 6, = 42 in 
(27) it is seen that the driving force always opposes the motion. 
Hence solvent circulation is possible only if the tubes have different 
structure resistance factors. 

Expression (15) may be written as 


CoC" ; 
A=B/t——e |; (30) 


and with (22) this becomes 
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fs 1 BL ce’ + ¢" 
A=B phe ha ee Aisa ee 
CSC poe ae |) @ 
or 
oy ed Fi ce’ + ¢" 
Aesth eG a ae ea 
(c’ 0") (—7—=-+--) —B—5 (82) 
Neglecting ¢ in comparison with fC, in (29) we have from (27) 
a L 
J Hes (A, 6,— As 62) , (33) 


and with (32) and (29) we find, to the first power of B in (32), 


a L (ou — Cc” 
HE ee af tanh UN eee ind 
76,(¢ 6") (62 — 4) (Ee | Ke + 02) : 
4B (6, ki — 02 fr) (84) 
Noting that the last term in the bracket is proportional to J this gives 
be+ec Caen ce J 
cca ae a 6) 
L , ” 

= fen +") (0, -+ 2) . (35) 


Since be << C, in (7) the left side is very nearly equal to J , and (6) 
then becomes 


KL E + (6+ 62) S| = (c’ —c”) (02 — 61) (36) 


The left side is approximately KL ©, and thus 


K= (6, — 61) (¢’ —c") . (37) 


104 
0 L{Cy 
The signs are such that the solvent flows in the same direction as the 
solute in the tube where the structure 1s more porous. The speed is 
proportional to 


D6é1cé—ec" 
ke RE a 38 
0M L EEL 
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ERRATA 


In the paper “The Mechanism of Cell Division” by N. Rashevsky, 
pp. 23-30, the following somewhat confusing notation should be 
changed: substitute J for T in last term of equation (1) ; in left hand 
side of equation (3); in first line from top on page 25; in 12th and 
19th line from top on page 27; in left hand side of equation (20). 

In formula on 10th line from bottom on page 24 the expression in 
braces should be multiplied by c. 


In the paper “On the Mechanics of Viscous Bodies and Elongation 
in Ellipsoidal Cells” by Gale Young, pp. 31-46, replace V by v on page 
31 last line; page 32, equation 4; page 33, lines 4 and 5. In equations 
(6) and (9) replace y by 7 in arguments of the cosine function, pages 
29.22. Also in equation (13). ” denotes direction of normal to sur- 
face. Change ¢ to ® in equation (32) and 6th line from bottom of 
page 37. Replace n by 7 on right side of (63). 


In the paper by Gale Young, “Blongation and Constriction in 
Cell Division,” pp. 75-91, in equation (36) page 85, replace (cos nz) 
by cos (nz). Change v to V in equation (48) page 86. 


In the paper by N. Rashevsky, “Mathematical Biophysics of 
Growth,” pp. 119-127, in equation (26) page 123, the coefficient of the 
log term should be 3/2. 
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